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ABSTRACT 



A procedure is developed which can be used to obtain the 
solution of a two-sided allocation problem with payoff of the 
form I = 5q K(x,y,t)dt if K is concave in x(t), i.e, K < 0, 

and convex in y(t) for all t in the interval [a,b], A 
second degree polynomial example is explained in detail to 
illustrate the procedure. The computer program and a sample 
output for this example constitute the Appendix. Next an 
exponential search problem is discussed. Then in conclusion 
a general application of the method is outlined. 
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1, Introduction., 

The purpose of this thesis is to investigate a method for 
obtaining the optimum solutions of two-sided allocation prob= 
lems. A two-sided allocation problem is one in which there 
are two sides allocating resources with opposing interests. 
Generally one side can optimize its position by minimizing 
some factor, such as 50 ^^"^ where each side has some control 

Oi 

of K, which the opposing side would like to maximize, or vice 
versa. 

The problem of maximizing (or minimizing) a functional 
K[y(t)] dt over a certain class of functions (y(t)] has 
been studied by use of various methods. Among these are the 
use of Lagrange multipliers[2, 9], techniques of calculus of 
variations[l , 8], and methods of control theory[3] such as the 
differential methods used in this paper. 

There appears to have been little application of these 
methods to the two-sided allocation problem where K is a 
fmction of two functions, x(t) and y(t). Although the theory 
of games has application to this type of problem if K(x,y,t) 
is a real valued f motion that is convex in y, i.e, K > 0, 

*/ J 

and concave in x[6], it does not provide a method for obtain- 
ing a solution. The subscript indicates partial differentia- 
tion with respect to the corresponding variable. 

Until the last five years there did not exist a method 
which could be used to obtain the solution of optimization 
problems subject to constraints on resources. However, differ- 
ential methods such as those used by Faulkner[3] have been 

% 
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developed recently which can solve such problems ^ A method 
for determining solutions by differential equations using an 
analog computer has been developed by Kose[4]o 

A version of the differential methods given by Faulkner 
[3] will be used in this thesis* These methods not only seem 
simpler to grasp and apply to the particular functionals con- 
sidered, but they provide a method for solving either maxi- 
mization, minimization, or minimization -maximization (saddle 
value) problems with a speed and simplicity that none of the 
others (with the possible exception of Kose[4]) seem to have. 
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2, Basic principles. 

Consider a functional of the form 

(1) I = K dt 

where: K = K(x,y,t) is a real valued function; x = x(t) is a 

fmction to be found which maximizes (1) subject to 

(2) X > 0 and 5^ x dt = 1: 

and y = y(t) is a function to be found which minimizes (1) 
subject to 

(3) y > 0 and 5^ y dt = 1, 

A function x, or y, which is piecewise continuous in [a,b] 
and satisfies (2), or (3), will be called admissible. If the 
total resources available to each side are not unity, they 
can easily be normalized to comply with (2) and (3), making 
the functions x(t) and y(t) density fimctions which indicate 
relative allocation of resources. 

It will be convenient to introduce a new functional, J. 
Let H = K + Ax +^y where A and ^ are constants (Lagrange 
multipliers) to be found. Let 

(4) j = 5^ H dt. 

We see that if x and y are admissible functions, then J dif- 

fers from I by the number A + /<. 

Now if y were a given fimction we could get a maximum 
value for J, 

J(xo,y) = maximum J(x,y) 

by maximizing the Hamiltonian. H, subject to the constraint 
X > 0, for each value of t on the interval [a,b]; however 
the fimction x(t) might not be admissible. Also, for given 
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x(t) we could get a minimum 



J(x^yo) = 



minimum 

y 



J(x,y) 



It should be pointed out that the function Xq above depends 
on both y(t) and the parameter A, and yo depends upon x(t) 
and yUo Observe that also 

(5) J(x,yo) < J(xo,yo) < J(xo,y) 

since by the way it was obtained x© maximizes J(x,y) for any 
given y, even y©, and yo minimizes J(x,y) for any given x, 
even Xq . 

If the functions Xq and yo are admissible, we will call 
them the min-max strategies (x’^^jy^) and we will call J(x*,y*) 
the min -max solution o Hence, to solve a min -max problem we 
want to determine, for some A and yU initially unknown, 

J(xo?yo) subject to the restrictions (2) and (3) on x(t) and 
y(t)o Any pair of admissible functions (x^,y’**’) which fur- 
nishes the min-max solution of (4) will also provide the min- 
max solution of (1) since 

(6) J = H dt = 5^(K + Ax +/<y)dt = I + (A +/<)» 

Let us assume that K has continuous second derivatives » 

Then notice that an absolute maximum (or minimum) which is 
not a stationary point can only occur at x = 0 (or y = 0) for 
admissible functions^ If an absolute maximum occurred for 
X = Xi >0, then any larger value of x, say x = Xj + £, would 
yield a larger value of H since H > 0 at Xj . A similar argu- 
ment applies to minimization over y» Therefore for a min-max 
solution to exist, the functional K(x,y,t) must not be mono- 
tone increasing in x or monotone decreasing in y» 
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If the second derivatives of K are continuous, then the 
follo\’/ing conditions exist. To get maximum H at a point on 
the interval [a,b], either there is a relative maximum for 
which 

(7) " ^x ^ ® 

or else, if the maximum occurs on the boundary then 

(8) X = Oo 

Similarly for minimum H, either there is a relative minimum 
for which 



(9) Hy = Ky +/< = 0 

or else, if the minimum occurs on the boundary then 

(10) y = 0, 



An indication of a boundary situation will be negative values 
of X or y as solutions of equations (7) and (9), 

A min-max problem may also be called a saddle value 
problem if the integrand is strictly concave in x (K < 0) 
and strictly convex in y (K > 0) on the interval [a,b], 

•X «/ 

According to Kuhn and Tucker [ 6] the conditions 

H^Cxcyo) < 0, H (xo,yo)xo =0, Xo > 0 

( 11 ) ^ ^ 

Hy(xo,yo) ^ 0, Hy(xo?yo)yo = o, yo > o 

are necessary in order that (xo,yo) provide a solution for 
the saddle value problem, which neatly summarizes the previous 
two paragraphs. 

If the function K is not concave in x and convex in y 
the possibility of multiple solutions of (7) and (9) arises. 

In this case the difficulty of computation increases consider- 
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ably. Each possible pair of roots (x,y) must be checked to 

determine the x which maximizes H and the y which minimizes 

H, This operation must be repeated for each value of the 

parameter t on the Interval [a,b]. 

Only the simple case where K(x,y,t) is strictly concave 

in X, K < 0, and strictly convex in y, K,_ > 0, over the 
XX yy 

domain [a,b] of the parameter t will be considered further. 
Differentials will be defined and a Newton-Raphson iteration 
set up to determine the desired solution. 

While we had no trouble with convergence, this is often 
a problem when using a Newton-Raphson iteration procedure. 
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3, Computational procedure. 

In this section we will set up a Newton -Raphson iteration 
for the solution. To start, let us guess values for the La- 
grange multipliers A and yU and compute a trial solution for 
x(t) and y(t) using equations (7), ( 8 ), (9) and (10), Usually 
it will not be an admissible solution since 

5^ X dt = X 1 and 5^ y dt = Y / 1, 

But we shall see that if the solution thus found is, by rare 
good luck, admissible, then we have the min-max solution. 

This follows directly from (5) since any trial solution is 
(xo^yo) for the particular values of A and yU chosen. 

Let us now set up a routine for correcting the values of 
(X,Y) obtaining a sequence approaching (1,1), Suppose we con- 
sider the given trial solution and neighboring values of A and 
say A + dA and /a. + c^. Every optimal solution must satisfy 

H„(xo,yo)xo = 0 and H (xo,yo)yo = 0 from equations ( 11 ), Gen- 
^ y 

erally only one of the terms in each product is zero at any 
point on the interval [a,b]. If these zero terms are con- 
strained to remain zero as we change x, y, A and by small 

amounts (fx, <5y, dA and (^, then x and y must always maximize 
and minimize, respectively, H since the requirements (11) re- 
main satisfied. We must change x, y, A and ^ in such a way, 
satisfying the above, that the values X and Y approach unity. 
In order to calculate 

dX = 5^ 6k dt and dY = 5^ 6y dt 

cL 3 . 

let us divide the interval [a,b] into four subsets of inter- 
vals To, Ti , Tg and T 3 , On the intervals in Tq the min-max 
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values of x and y are obtained from (7) and (9), i.e. x > 0 
and y > 0; on those in Tj we use (7) and (10), i.e. x > 0 and 

y = 0; on those in Tg we use (8) and (9), i.e. x = 0 and 

y > 0; and on those in T3 we use (8) and (10), i.e. x = 0 

and y = 0, The resulting variations c5x(t) and ^y(t) will be 

given in the intervals of T© by 



( 12 ) 



0 = (Ty + dA 

0 = <5h = H rfx . H <fy + d//. 



y yx -- -yy 
In the intervals of they are given by 

0 = <^H = H (fx + dX 

(13) ^ ^ 

0 = dy. 

since y = 0 on these intervals. In the intervals of Tg they 
are given by 

0 = Sx 



(14) 



0 = jHy = 



H. 



yy 



dy + d/t 



(16) 



since x = 0. Finally in the intervals of T3 
0 = 6x 

(15) 

0 = dy. 

Hence we have equations of the form 

bii(t) Sx + bi 2 (t) cfy + dA = 0 
b2 i(t) d"x + b 2 g(t) <5y + d/i= 0 

which can be solved for 

dx(t) = aii(t) dA + ai 2 (t) d/t 
cTy(t) = a 2 i(t) dA + a 2 g(t) dyU 
where the are obtained from (12), (13), (14) or (15), 

depending on the interval. 

We are now ready to calculate the total effects, dX and 



(17) 
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dYj of the small changes dx{t) and c5^y(t)„ Integration over 
the entire interval [a,b], consisting of the four subsets Tqj 
T j, Tg and T 3 , yields 

dX = All d^ + Ajg 
dY = Agi dA + Agg dyu 



(18) 

where 



dX = 5^ cTx(t) dt, 

3 . 



dY = (^’y(t) dt. 



Now let us treat the differentials in (18) as differences to 
get 



Ai 1 aA + Ai 2 A>^ — AX “ 1 — X 

(19) 

Ag 1 A A Agg — AY — 1 “ Y„ 

These equations (19) can now be solved for AA and /yt; thus we 
determine new values of A andy^ which tend to drive the solu- 
tion toward admissibility. 

Let us repeat this. We compute another trial solution 
using the new values of A end Repeat this sequence until 
the "errors" dX and dY become sufficiently small and the 
solution is admissible, within the required accuracy. 

It is convenient to define an error term which can be 
used to check the convergence toward admissibility. Let the 
error, E, be defined as 

(20) E = (1 - X)^ + (1 - Y)^ = AX^ + aY^, 

AS the error goes to zero the solutions approach admissibility 
since both X and Y must be approaching unity and the differ- 
ences, AX and AY, must approach zero, 

important feature of this differential method is that 
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each solution is an optimal (maximizing and minimizing) solu- 
tion for the particular values of 'X and , or endpoints X and 
Y, considered o The successive iterations drive Xq and Yq to 
admissibility j thus yielding the min-max (or saddle value) 
solution v/hich satisfies the given constraints (2) and (3)o 



10 = 



II 



4. Second degree polynomial examples. 

In order to test the method a computer program was 
written, using FORTRAN language, and run on a CDC 1604 com- 
puter. The form of this program and a sample solution, which 
will be explained in detail, appear in the Appendix. The 
program was set up first for an integrand of the type 
(21) K(x,y,t) = (y + At + B)^ + xyt - (x + Ct + D)^. 

The parameters A, B, C and D were assigned various sets of 
integer values to test the general applicability of the 
method. In every case the interval [a,b] was chosen to be 
[ 0 , 1 ]. 

The choice of initial values of the Lagrange multipliers, 
A and^, is at times a minor problem in getting started. If 
the initial values yield only negative values of x(t) and y(t) 
over the Interval [a,b], as determined by equations 
(7) H = + A = 0 and 

(9) Hy = Ky = 0, 

then the procedure has no way of determining appropriate 
changes, AA and of the Lagrange multipliers for a new 
iteration since each of the A . . ' s in equation 

dX — Aj ^ dA ^12 
(18) ^ ' 
dY = A 21 dA + A 22 

will be zero. The solution of equations (7) and (9) at any 
point on the interval [a,b] with the requirement that both x 
and y be positive at that point should Indicate appropriate 
initial values of A and /A such that convergence and a solution 
will occur. 
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Now let us go through the method in detail for a specific 

(22) K(x,y^t) = (y + 2t)^ + xyt = (x + 3t)^„ 

The Hamiltonian, H, is 

(23) H = (y + 2t)^ + xyt - (x + 3t)^ + ^x + /^y 
and its second partial derivatives are 



(24) 



The first partial derivatives of H are 

E = -2x + ty => 6t + A 
(25) ^ 

Hy = tx + 2y + 4t + 

vVe can see that K is strictly concave in x and strictly convex 



H 


ct 

II 






xy 


xy 






H ^ = 


“2 = L 


< 


0 


XX 


XX 








2 = 


> 


Oc 


yy 


yy 







in y from (24) 
x( t) 



Upon setting = H^. = 0 we can solve for 




(26) 



y(t) 



6t‘ 



8t 
2 

t^ + 4 



At ^ 2yU 



Recall that the maximizing function x = Xq and the minimizing 

function y = yo satisfy 

Xq = x(t) if x(t) >0 

- 0 if x(t) < 0 

yo = y(t) if y(t) > 0 

= 0 if y(t) < 0 

when x(t) and y(t) satisfy equations (26) « 

To determine appropriate initial values of A and yu in 

this problem, evaluate (26) at a particular value of t, say 

t = Oo Then x(0) = 0,6A and y(0) = and any values of 

12 „ 



“X > 0 and yU < 0 are satisfactory for an initial guess. If we 
take ^ = 10 and /d ~ -10, then x(0) =5 and y(0) =5, 

Let us arbitrarily divide the interval [0,l] into 100 
equal parts of length 0,01 and calculate x(t) and y(t) by (26) 
for t = 0 , OoOl, 0 , 02 , » 0 ® , 0.39, 1,00, Meanv/hile we 
should determine which of the four subsets of intervals, Tq , 

Tj , Tg or T 3 , we are in for each value of t in order to cal- 
culate the a. .“s which are necessary for the determination of 
aA and ^ if the solution is not admissible. The a . . ' s are 
given: in Tq (x > 0 , y > 0 ) by 



1 ( t) “ 



K 






(27) 






^xxl^yy ^xy 
K. 



2 ’ 






in Ti (x > 0 , y = 0 ) by 



yy ° ^xy 



2 ’ 



2 ( ^ ) “ 






K 






^xx^yy " ^xy 



2 ’ 



K 



XX 



1 ( t) - 



(28) 



K ^ 

XX 



aig(t) - agj(t) - agg(t) 
in Tg (x = 0 , y > 0 ) by 

3-ii(t) = ajg(t) = agi(t) 



= 0 ; 



= 0 , 



(29) 



- K I 

yy 



K K - K 
XX yy xy 



2 ’ 



and in T 3 (x = 0 , y = 0 ) by 

(30) aij^(t) — ajg(t) = S 2 i(^) ~ 3-2 2 ^^) ” 

For this particular example we find that Tq = [O, 0,82], 
is empty, Tg = (0,82, 1 , 0 ] and T 3 is empty. 

Let us integrate Xo(t) and yo(t) using 
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a trapezoidal or. 



if more accuracy is desired, Runge-Kutta method to determine 
if they are admissible, i.e. they integrate to unity. They 
are not admissible since 



Xo dt = 4.1225* •• and 5 q Vo = 3,0644* 

For each value of t on [0,l], a change in x, y, A and^ ac- 



cording to 




^rx = 

(IV) 

oy = 


ail dA + ai2 d/i 

asi dA + a22 d/.i 



will retain the maximizing and minimizing properties of X£,and 
y^. The total changes over the interval [0,l] are 
dX = All dA + Ai 2 d/t 

(18) 

dY = A21 dA + A22 d/i 



where 




dX = 


(^x dt , dY = si (fy dt , 


A. . 


= dt. 



For first order approximations treat the differentials as 
differences and set 

AX = 1 - X = All AA + Ai2 = -3.1225* •• 

AY = 1 - Y = A 21 AA + A 22 - -2.0644* •*. 



Solve for 




aA = 


-5.3532* • • 


A/t = 


5 , 7408 • • • . 



Repeat the calculations of x and y by (26) v/ith 
A = 4. 6468* • • 
yU =-4. 2592* * • . 

Continue the routine and repeat until AX and aY are acceptably 
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small. 



Other polynomials solved by this procedure were; 

(31) K(x,y,t) = [1 + (t - .5)^]y^ + xyt - x^, 

(32) K(x,y,t) = y^sin(t) - xycos(t) - x^. 

The last one was solved for the interval [0,l] and also the 
interval [O,ir]o Care must be taken to ensure that the inte= 
grand, K, maintains its strictly concave, strictly convex 
properties over the entire interval. The interval [0,2ir] 
cannot be used for (32) because K = 2sin(t)j which becomes 
negative on the interval [ir,2Tr]. 
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5. Exponential search exampleo 

The particular functional which kindled interest in the 
two-sided allocation problem was the integral 

5^ x[l = e“^^] dt 

which represents the probability of detecting an object dis- 
tributed as X = x(t) over the interval [a^b] by use of a 
search effort distributed as y = y(t) over the same interval 
if the law of detection has the exponential- saturation form 
as in random search [?9 9], The function g = g(t) is a measure 
of the ease with which a unit of search conducted near t 
detects a target there and is considered known., For "stand- 
ard" conditions^ g may be considered unity„ 

To maximize this integral with respect to y it is neces- 
sary to minimize over y 

(33) I = 5^ X e“^y dt„ 

oL 

This problem was posed and a method of solution given in 
Search and Screening ., OoEoG„ Report No. 56[9] and in a paper 
written by Bernard 0. Koopman[7] who edited the above report. 
However j the maximization over x of 

min 5^ X e“®^ dt 

y 

was not considered. 

Notice that this problem does not meet the requirement 
of strict concavity in x since K = 0. However the differ- 
ential method produced a min-max solution and min-max strate- 
gies with no difficulty. Minor modifications to the FORTRAN 
program listed in the Appendix were made to fit it to this 



problem. The function g was assumed to be of the form 
g = A + Bt + Ct^ + Dt^ 

and required to be strictly greater than zero. The parameter 
A was given the values 0,01, 0,1, 1, 3 and the other parame= 
ters were either all zero or else one was set equal to ±0,9 
until 17 different problems were solved. The smallest and 
largest values of g were 0,01 and 3,0, The number of itera- 
tions required to drive the error less than 10” was never 
more than six. The program compiled in 47 seconds and required 
about 10 seconds to complete a single solution. 

It became apparent that this particular problem could be 
solved more easily than by the differential method used, Con» 
sider the Hamiltonian 

H = X + Ax + /<y. 

To minimize H over y(t) we solve 

H = »gx = 0, 

V 

We may solve for y, obtaining 

y(t) = ~ log if g(t)x(t) >// 

(34) S e 

= 0 if g(t)x(t) <yUo 

This is the solution obtained by Koopman[7, 9] for any fixed 
function x( t) „ We must also maximize H over x(t) and hence 
we have a second equation 

+ A = 0, 

We must solve these simultaneously for x and y„ We require 



then 



or g(t)y(t) = constant „ From (34) we see that gx must also 
be constant if gy is to be constant „ Due to the requirements 
for admissibility it follows that 

x(t) = y(t) = 

where T is a constant to be determined „ If the constant T is 
assigned any reasonable positive number, one iteration will 
be sufficient to determine the min-max strategies (x^,y*) 
since the initial solution (xo,yo) can be made admissible by 
normaliZ-ing ito For example if 

5^ Xo dt = X = 5^ yo dt 

then 

Xaitl ^ _ ^b ynlt) 

3 , A 3 A 

and 

X^'ct) = - y'^(t)o 

All that remains is to show that the solution attained 
is truly a min=-max solution in that x*^ and y* are min»max 
strategies „ Consider the payoff 

I(x%y") = 5a x"e“^y"dt = 5^ - A xMt = » . 

Observe that if gy^ remains constant over the interval [a,b] 
then for all admissible x we have 

Kxjy'^) = dt = = I(x'',y^)„ 

Since the integrand, K, is strictly convex in y (K >0), it 

V •/ 

is obvious that for all admissible y(t) 

Kx'^jy'^) < I(x*^,y) „ 

If X - Xi(t) / x’^(t), then there exists a yi^(t) such 

18 o 



that 



I(xi,yi^) < I(x%y") 

and X x**(t) obviously does not maximize I, In conclusion 
Kxijyi’*) < I(xi,y^) = I(x’^,y*') < I(x^,y) 
and (x’^^y^) is a roin-max strategyo 

As a rough check to ensure that x^ was a maximizing stra 
tegy, equation (34) was used to calculate the minimizing stra 
tegies y for predetermined non=optimal strategies 
x(t) = 1, 

x(t) = 0o06 + 1,91, 
x(t) = lo95 “ 1.9t, 

letting g(t) vary as in the min-max calculations previously „ 
In every case the min-max solution payoff was greater than 
the minimizing solution payoff when x was fixed not equal to 
x^(t)« 
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6. Conclusions o 



The differential methods discussed in this thesis can be 
used to determine 



mgx min 5^ K(x,y,t) dt 



subject to 



X > 0, 



(35) 



5a X dt = X 
y > 0, 5^ y dt = Y 



where X and Y are given numbers and the integrand has certain 
propertieSo It must not be monotone increasing in x or mono- 
tone decreasing in y^ for x and y large, over the interval 
[a,b]« The integrand, K, is assumed to be finite and real 
valued although it need not be continuous as long as the 
regions of discontinuity are known* Determination of maxima 
and minima is accomplished by comparison of the value of K at 
stationary extrema with its value at the discontinuities* 

This thesis specifically covers the cases when K is twice 
differentiable, strictly concave in x (K < 0), and strictly 

convex in y (K > 0) for non=negative piecewise continuous 
yy 

X and y over the domain [a,b] of the parameter t. For this 
type of function, K, a min=max solution, Kx'^(t), y^(t)), and 
a min=-max strategy (x*^(t), y^(t)) can be determined which 
satisfy the constraints (35)* 
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appendix 



First the symbols and terms "used in the FORTRAN program 
will be correlated with those introduced in the text of the 
thesis; then a brief explanation of the program, the program 
itself, and a sample output will follow. 



t^ORTRAN 

symbols 


Text 


Explanation 


T 


t 


the independent variable defined 
the interval [0,l] 


XM 


xo 


that X which maximizes H where 
H = K + X X + JAY 


YM 


7o 


that y which minimizes H 


FUNF 


K 


I = 5^ K dt where K = K(x,y,t) 


FXYF 


^xy 


a^K/exsy 


FXXF 


^xx 


a^K/ax® 


FYYF 


V 


a^K/ay^ 


A(l,l) 


1 ( t) 


from equation (17) 


AT0T(1,2) 


at 


from equation (18) 


U 


A 


Lagrange multiplier 


V 


/A 


Lagrange multiplier 


XTOT 


X 


5^X0 dt 


YTOT 


Y 


5Jyo dt 


T(I) 


t 


a particular value of t on [0,l] 


X 


x( t) 


a solution of equation (26) 


Y 


y( t) 


a solution of equation (26) 


DX 




a small change of Xo(t) 


DY 




a small change of yo(t) 


DU 


dA 


a small change of /\ 
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FORTRAN 

symbols 


Text 


Explanation 


DV 




a small change of /A 


DIFFX 


dX 


1 - 5 ^Xo dt 


DIFFY 


dY 


1 " 5oYo dt 


PAYOFF 


I(x,y) 


5 ^K dt 
•^0 



Items set off by a brace of CCCCCC's will generally 
change when the Integrand, K, is changed. 

After the maximum number of Iterations, division of the 
interval [a,b], acceptable error and Initial Lagrange multi- 
plier values are determined, we commence the iterations. The 
interval [a,b] has been chosen to be [0,l]. The variable T(I) 
is defined In such a way that it takes on values from zero to 
one as the index, I, goes from 1 to 101. We now solve the 
first derivative equations (26) for X and Y, and determine 
which set of subintervals, T 3 , T 2 , Tj or Tq , we are in so that 
XM, A(l,l), A(l,2) , A(2,l) and A(2,2) may be determined. 

Since we are using a trapezoidal Integration routine, we 
sum XM, YM, A(l,l), A(l,2), A(2,l) and A(2,2) as we calculate 
them and then multiply each total by the width of one inter- 
val, At, to determine the integrals. After the error is 
determined, pertinent values are printed and then if the solu- 
tion is not admissible, within the acceptable error, changes 
of the Lagrange multipliers are computed using the subroutine 
at the end of the program. The iterations are repeated until 
the error is acceptable and the solution XM, YM is admissible. 
The payoff is now calculated using again the trapezoidal 
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integration method „ The min-max solution (PAYOFF) and min- 
max strategies (XM, YM) are now printed. 

Statements 160, 140 and 135 are used only if poor ini=- 
tial values of the Lagrange multipliers are chosen such that 
the matrix ATOT is singular (i,e, its determinant is zero). 
The printout is for information only. 

The numbers at the end of the program are sample values 
of AA, B, C and D to be read into the program as data. 
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I 



PROGRAM MINMAX 

DIMENSION TI900), XM(900), YM ( 900 ) , . ERROR ( 50 I , 
1 A(2,2), ATOT(2,2> 

C 

C INSERT THE FUNCTION TO BE INTEGRATED AND ITS 

C SECOND PARTIAL DERIVATIVES 

CC 

CCCC 

CCCCCC 



FUNF(X,Y,T> = (Y+AA»T+B)»«2+X*Y»T-(X+C*T+D)**2 
FXYF(X,Y,T) = T 
FXXF(X,Y,T) = -2. 

FYYF(X,Y,T) = 2. 

CCCCCC 

CCCC 

CC 

C MAXIMUM NUMBER OF ITERATIONS = NR. .MUST AGREE WITH 

C ARRAY SIZE OF ERROR 

NR = 50 
C 

C NUMBER OF INTERVALS OVER WHICH TO EVALUATE 

C AND INTEGRATE 

NIPLUSI = 101 

C I 

C EPSILON IS THE ACCEPTABLE ERROR 

EPSILON = . lE-15 
C 

WRITE OUTPUT TAPE 6, 200 

200 FORMAT! 1H1,////////31X, 18HS0LUTI0N OF MINMAX////) 
CALL TIM£(6) 



C 

CC 

CCCC 

CCCCCC 

2 

8 

210 

1 

2 

CCCCCC 
CCCC 
CC 



READ INPUT TAPE 5, 8,AA, B, C, D 
FORMATIUFIO.S ) 

WRITE OUTPUT TAPE 6, 210, AA, B, C, D 
FORMAT! 16HK(X,Y,T) = !Y + ,F8.5,5H*T + ,F8.5, 

19H)**2 + X*Y*T ~ !X +,F8.6f 5H *T + f-F8.5, 
UH)**2///) 
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1 



C 
C 
C 
C 

CC 

CCCC 
CCCCCC 



WRITE OUTPUT TAPE 6, 220 

FORMAT! 3HRUN, 1 OX , 1 HU, 1 6X , 1 HV , 15X,4HXT0T, 13X,4HYT0T, 
14X,5HERR0R// ) 

SET THE INITIAL VALUES OF THE LAGRANGE MULTIPLIERS 



U = 10. 

V = -10. 



CCCCCC 

CCCC 

CC 

c 

FNI = NIPLUSI - 1 
C 

C COMMENCE ITERATIONS 

C 

DO 100 L = 1 , NR 
XTOT = 0. 

YTOT = 0. 

DO 5 M » 1.2 
DO 5 N = 1 ,2 
5 ATOT!M,N). = 0. 

C 

C COMMENCE CALCULATIONS OVER THE INTERVAL !0,1) 

C 
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i 



c 
c 

cc 

cccc 

cccccc 



DO 99 I = 1 , NIPLUS 1 
FI = I - 1 
T( I ) = FI/FNI 

SOLVE FIRST DERIVATIVE EQUATIONS FOR X,Y 



1 

CCCCCC 

cccc 
cc 
c 



DENOM = T( I ) **2 + 4. 

Y = (T( I )*(2.*C*T( I ) - 4.*AA + 2.*D - U) - 4.*B 
1 - 2.*V)/DENOM 

X = (-T(I )*(2.*AA*T( I ) + 2.»B + 4.«C + V) - 4.*D 
+ 2.*U)/DENOM 



C 

C 

C 



C 

C 

C 



C 

c 

c 



c 

c 

c 



c 

c 

c 

c 

c 

c 

c 

c 



IF (X) 10,20,20 
10 IF (Y) 15, 25, 25 
20 IF (Y) 35, 45, 45 

BOTH X AND Y NEGATIVE 

15 XM( I ) = 0. 

YM{ I ) = 0. 

A( 1 , 1 ) = 0. 

A( 1 ,2) = 0. 

A(2 , 1 ) = 0. 

A (2, 2) = 0. 

GO TO 50 

X NEGATIVE AND y’ NON-NEGAT I VE 

25 XM( I ) = 0. 

YM( I ) = Y 
A( 1 , 1 ) = 0. 

Ad ,2 ) = 0. 

A(2, 1 ) = 0. 

A(2,2) = -1 ./FYYF(XM( I ) , YM(I), T(U) 

GO TO 50 

X NON-NEGATIVE AND Y NEGATIVE 

35 XM( I ) = X 
YM( I ) = 0. 

A(l,n = -1 ./FXXF(XM( I ) , YMin, T( I) ) * 

A( 1 ,2 )= 0. 

A(2, n = 0. 

A(2,2) = 0. 

GO TO 50 

BOTH X AND Y POSITIVE 
45 XM(I) = X 

YM( I ) = Y ^ 

DET = FXXF(XM( I ) ,YM( n , T< I ) )*FYYF( XM ( I ) , YM( n ,Tl 
1 -FXYF(XMd), YM(I), T(in»*2 ^ 

A(1,l) =-FYYF(XM(n, YM( I ) , • T{ I n/DET 
A (1,2) = FXYF(XM(I), YM ( I ) , T(I))/OET 
A(2, 1 ) = A( 1 ,2) 

A(2,2) = -FXXF(XM( I ), YM(I), Td))/OET 



50 

52 



DX = Ad , 1 )»DU + A( 1,2)*0V 
,DY = A(2,1)*DU + A(2,2)»DV 

TRAPEZOIDAL INTEGRATION ROUTINE 
IFd-1 )52,52,54 

XTOT=XTOT+XM( n/2. ». 

YTOT=YTOT+YM( I )/2. 
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ooo oooooooo ooo oooo 



DO 53 M=l,2 
00 53 N=1,2 

53 AT0T(M,N)=AT0T(M,N)+A(M,N)/2. 
GO TO 99 

54 IF{ I-NIPLUSl )55»52,52 

55 XT0T=XT0T+XM( n 
YT0T=YT0T+YM( I ) 

DO 56 M=1,2 

00 56 N = 1,2 

56 AT0T(M,N)=AT0T(M,N)+A(M,N) 

GO TO 99 

99 CONTINUE 

XTOT = XTOT/FNI 
YTOT = YTOT/FNI 
DO 60 M = 1,2 
00 60 N= 1,2 

60 AT0T(M,N) = ATOT ( M, N 1 /FN I 



DETERMINE ERROR 

OIFFX = 1 . - XTOT ' ; 

DIFFY = 1 . - YTOT ' . 

ERROR(L) = 01FFX**2 + DIFFY**2 

WRITE OUTPUT TAPE 6, 2 30 ,L ,U ,V , XTOT. YTOT , ERROR ( L ) 

230 F0RMAT(I2,3X.F14.8,3X,FlU.6,5X,F12,e ,5X,F12.8,5X, 

1 E 1 4 • S J ’ 

IF(ERR0R<L) - EPSILON) 105, 95, 95 

CALL SUBROUTINE TO SOLVE MATRIX FOR. OU AND DV 

95 CALL SOLVE! AT0T( 1 , 1 ) , ATOTI 1 ,2) ,ATOTt 2, 1 1 ,ATOT(2,2 ) , , • 
1 DIFFX, DIFFY, DU, DV, IND) 

GO TO (97,150), IND 
97 U = U + OU 
100 V = V + OV . 



END OF A RUN 

REPEAT CALCULATIONS USING THE NEW VALUES 



OF U AND V 



BEST U, V, XM(I), YM(I) ARE NOW STORED FOR EACH T { ,I ) 
INTEGRATE FUNF(XM(I), YM ( I ) , T(I)) 

105 PAYOFF=0. 

DO 70 I=1,NIPLUS1 
IF{ I-l 162,62,64 

62 PAYOFF = PAYOFF + FUNF(XM(I), YM(I), T(I))/2. 

GO TO 70 

64 IF( I-NIPLUSl 165,62,62 

65 PAYOFF = PAYOFF + FUNF(XM(I), YM(I), T(I)) 

70 CONTINUE 

PAYOFF=PAYOFF/FNI 

END OF COMPUTATION 

WRITE OUTPUT TAPE 6, 110, PAYOFF 
no F0RMAT(//9HPAY0FF = ,E14.8//1 

WRITE OUTPUT TAPE 6,130 ( T ( I ) , XM ( I ), YM ( I) , 

1 I = 1,NIPLUS1) 

130 F0RMAT(2TX,lHT,18X,2HXM,20X,2HYM,//( 1 7X , F8. 5 , 1 OX, 

1 F12.9, 10X,F12.9) ) 

GO TO 160 

150 WRITE OUTPUT TAPE 6,140 

140 F0RMAT(5X,31HDETERMINANT OF THE MATRIX ATOT , 

1 IIHEQUALS ZERO///) I 

WRITE OUTPUT TAPE 6 , 1 35, XTOT , YTOT, U, V, ATOT( 1 , 1 I ♦ ' 

1 ATOT( 1 ,2 ) , ATOT( 2* 1 > tATpT(2,2 ) 

135 F0RMAT(7HXT0T « , E 1 2. 6 , 1 OX, 7HYT0T » , E 1 2 . 6/ 3X, 4HU - , 

1 E12.6, 13X,4HV = , E12.6//1X.6HA1 1 = ,E12.6, 

2 11X,6HA12 = ,E12.6/1X,6HA21 = ,E12.6,11X, 

3 6HA22 = ,E12.6/) 
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unuluiro 



160 CALL TIME(6) 

WRITE OUTPUT TAPE 6t 240 
240 FORMAT(lHl) 

GO TO 2 

C READ DATA FOR ANOTHER PROBLEM 

END 
C 
C 

SUBROUTINE SOLVE! A1 1 , A 1 2 , A2 1 ♦ A22 , B 1, B2, X 1 ,X2 , IND ) 
C 

C SOLVES THE MATRIX AX=B FOR X ' 

C 

D=A1 1*A22-A12«A21 
I F ( D ) 1 , 1 2 . 1 
1 Xl=A22/D*Bl-A12/D*B2 
X2=-A21/D*B1+A11/D*B2 
IND=1 
RETURN 
12 IND=2 . 

X1=0. 

X2 = 0. 

RETURN 

END 

END 




3. 

4. 
3. 
3. 



0 . 

1 . 

2 . 

2 . 
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SOLUTION OF MINMAX 
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